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A new next-to-leading order Monte Carlo program for the calculation of fully differ- 
ential heavy quark cross sections in electroproduction is described. A comparison 
between the theoretical predictions and the latest charm production data from HI 
and ZEUS at HERA is presented. 



1 Introduction 

Heavy quark production in deeply inelastic electron proton scattering (DIS) 
has long been recognized as a candidate process for extracting the gluon dis- 
tribution in the proton iecause the process is almost exclusively dominated 
by photon-gluon fusion EJ. The cross section for heavy quark production in 
DIS is given in terms of the heavy quark structure functions F2{x, , m?) and 
Fl{x,Q^ ,Tn^). A calculation of these in next-to-leading order (NLO) QCD 
has been available for several yearsH. However,^nly recently have they been 
calculated in a fully differential form at NLO B. The relation between the 
cross section and the structure functions, in terms of the usual DIS kinematic 
variables x, y and Q^, is as follows: 



dxdy Q* 

The electromagnetic coupling is a = e^/An, ^/S — 300 GeV at HERA, and m 
is the heavy quark mass. 

2 Method 



The r£sults presented herein are based on the fully differential structure func- 
tions u and Eq. |^. The computation was performed using the subtraction 
method. Therefore, all phase space integrations may be evaluated using stan- 
dard Monte Carlo techniques to produce histograms for fully differential, semi- 
inclusive, or inclusive quantities related to any of the outgoing particles. Ex- 
perimental cuts can also be imposed. When considering the production of 
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Figure 1: Charm structure function F2{x,Qni,m^) as a function of x for several values. 
Data is from DIS D" and D* production Q. Theory is NLO QCD as described in text. 

charmed mescms, fragmentation effects are modeled using the formahsm of 
Peterson et alu. 

3 Results 

Shown in Figure 1 is the charm contribution F2{x, Q^, m^) to the proton struc- 
ture function as derived from inclusive D*^{2Q1Q) and Z?°(1864) production 
by HI □ (circles) as a function of x at several values. The EMC El data 
are also shown (boxes). The solid lines are NLO QCD0 using GRVB piarton 
distribution functions. The dotted lines are NLO QCD using CTEQ 13. In 
both cases /i = /x/ = /i^ = \JQ'^ + 4to^. The sensitivity of the NLO calcula- 
tion with respect to variations in the renormalization (factorization) scale is 
greatly reduced relative to the LO calculation amounting to ±10% at small 
X when the scale is varied between 2m.c and S-^/Q^ 4- 4to^. At small the 
largest uncertainty is from the charm quark mass. A ±10% variation about 
the nominal value of 1.5 GeV used here results in a variation of ±15% at small 
X and = 12 GeV which decreases to ±8% at = 45 GeV. 

In Figure 2 D* production cross sections for the kinematic range 5 GeV^ < 
g2 < 100 GcV^ and < y < 0.7 with 1.3 GeV < pf * < 9 GeV and \ri°' \ < 1.5 
are shown. The data, shown as a function of (a) pj^ , (b) -q^ , (c) W, and 
(d) Q^, are from a recent ZEUS analysis □. The various lines are NLO JDCDH 
using GRVH distributions supplemented with Peterson fragmentation □. The 
dashed line is for fi = 2mc, rric = 1.35 GeV and e = 0.035 and the solid line is 
fi = 2^yQ^ + Ami, = 1.65 GeV and e = 0.06. The agreement between the 



Figure 2: Differential cross sections compared to ZEUSQ DIS D* data. 

data and the theory is quite reasonable at this early stage of analysis. 
4 Conclusion 

The differential distributions presented here required placing cuts on several 
different variables. The flexibility of the Monte Carlo approach has enabled 
us to calculate all of these simultaneously. With the analysis of the 1995 data 
still to come and the accumulation of more data in the future we look forward 
to incorporating charm production in deeply inelastic scattering into a global 
fit to constrain the gluon distribution at low x. 
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